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1 Introduction 



The study of the current through a system in contact with two reservoirs at 
unequal chemical potentials or at unequal temperatures is one of the most 
studied aspects of the theory of non-equilibrium systems [1, 2]. 

For the last decade, there has been an increasing interest in the study 
of the fluctuations of the current of quantum particles (fermions) through a 
disordered wire. It is now well established that the quantum statistics of the 
particles determines the distribution of the fluctuations of the current, and 
that in the metallic regime [3, 4], this distribution [5, 6, 7] is universal. More 
recent works have shown that the main property of the quantum nature of 
the particles which was responsible for these universal fluctuations is the 
Pauli exclusion principle [8, 9, 10, 11, 12]. 

Here we consider the symmetric simple exclusion process SSEP [13, 14, 
15, 16, 17] which is a stochastic model of classical particles with hard core 
interactions (and without inertia) which diffuse on a finite chain with open 
boundary conditions. The chain is in contact at its two ends with two 
reservoirs of particles at unequal densities [18, 19, 20, 21, 22]. The com- 
bined effects of the stochastic injection and removal of particles at the two 
boundaries and of the diffusive nature of the hard core particles produce a 
fluctuating current. We calculate the first four cumulants of the integrated 
current. Based on our results for these four cumulants, we give a conjecture 
for all the higher cumulants and for the whole distribution of the current 
fluctuations. 

The fluctuations of the current in exclusion processes is also a subject 
with a long history [23, 24, 25, 26, 27]. Most of the known results obtained so 
far concern infinite geometries [23, 24, 27] or systems with periodic boundary 
conditions [26, 28, 29, 30] (see [31] for the variance of the integrated current 
of the asymmetric simple exclusion process with open boundaries). 

Our paper is organised as follows: in section 2 we define the model and 
we summarize our results. In section 3, we show how the first two cumulants 
can be calculated from the steady state properties. In section 4, we write 
a hierarchy (see also Appendix C) for the correlation functions on which 
our approach is based. In sections 5 and 6 we solve this hierarchy, in a 
low density expansion, where at each order the hierarchy can be truncated. 
Appendix A gives a derivation of the Gallavotti-Cohen relation [32, 33] for 
the SSEP with open boundaries. Appendix B points out the analogy with 
multi-particle ruin problems. 
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2 Definition of the model and main results 

2.1 The symmetric exclusion process with open boundaries 

In the one dimensional symmetric simple exclusion process, each site i (with 
1 < i < A^) of a one dimensional lattice of N sites is either occupied by 
a single particle or empty. A configuration C at time t is therefore fully 
determined by N binary variables Ti{t), the occupation numbers of the N 
sites {Ti{t) = 1 if site i is occupied and Ti{t) = if site i is empty). In the 
bulk, each particle independently attempts to jump to its right neighboring 
site, and to its left neighboring site, in each case at rate 1. It succeeds if 
the target site is empty; otherwise nothing happens (this means that during 
time t and time t + dt with < dt 1, a particle at site i jumps to site i—1 
with probability (1 — Ti-i)dt, to site i + 1 with probability (1 — rj_|-i)dt and 
does not move with probability 1 — (2 — Tj-i — Ti+i)dt). At the left boundary 
particles are injected at site 1 at rate a and removed from site 1 at rate 7. 
Similarly at the right boundary, particles are removed from site N at rate 
f3 and injected at site A'" at rate S. 

For general values of a,/3,7, 5, a current of particles flows through the 
system and we want to study the fluctuations of this current. To do so, 
we denote by Q{t) the number of particles which have moved from the left 
reservoir into the system during time t (so Q{t) is the number of particles 
which have jumped into the system at site 1 minus the number of particles 
which have left the system from site 1). We want to calculate the distribution 
of the total charge Q{t) during a long time t. 

For finite the system has 2^ internal configurations C (each site can 
be either occupied by a particle or empty). Let Pt{C) be the probability 
of finding the system in configuration C at time t. As the dynamics is a 
Markov process, the evolution of the probability Pt{C) of finding the system 
in configuration C at time t can be written as 



where we have decomposed the Markov matrix into three parts, depending 
on whether when the system jumps from configuration C' to configuration C, 
Q{t) increases by 1,0 or —1. (the matrix Wq contains all the diagonal terms 
which are all negative as well as all the non-diagonal elements corresponding 
to moves which do not take place at the left boundary, i.e. do not change 
Q{t)). One way to determine the distribution of Q{t) is to calculate its 




(2.1) 
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generating function (z'^^*)). 

If we define Pt{C, Q) the probability that the system is in configuration 
C at time t and that Q{t) = Q, one has 



= Y^W,{C,C')Pt{C',Q- 1) + Wo{C,C')Pt{C',Q) 
at ^, 

+W^i{C,C')Pt{C',Q+l) 
Then the generating functions Tt{C, z) defined by 



satisfy 

dVt{C,z) 
dt 



E 

c 



Vt{C,z)= Pt{C,Q)z^ 

Q=—oo 



z PFi(C,C') + W^o(C,C') + -T^-i(C,C'] 

z 



If we introduce the matrix defined by 



1, 



M^(C,C') = z + VFo(C,C') + -Ty_i(C,C') 

it is clear from (2.4) that in the long time limit 

c 



(2.2) 



(2.3) 



Vt{C\z) (2.4) 



(2.5) 



(2.6) 



where ii is the largest eigenvalue of the matrix ■ So this largest eigenvalue 
\i fully determines the distribution of (5(i) in the long time limit [26]. 



2.2 Symmetries of // 

In principle, \x depends on six parameters: the input rates a, 7, 5 at the 
two boundaries, the fugacity z and the number of sites N . There are three 
symmetries in the system that leave [i unchanged: 

1. The left-right symmetry: if we exchange the roles of a, 7 and (5, /?, 
this has the effect of exchanging the roles of the left and of the right 
boundaries, and so the statistical properties of Q{t) are replaced by 
those of —Q{t). Therefore /x should satisfy 

li{a, 7, 6, (3, z, N) = /x(5, /3, a,j,-,N) . (2.7) 



4 



2. The particle-hole symmetry: instead of counting the number of par- 
ticles Q{t) entering at the left boundary, one can as well count the 
number —Q{t) of holes entering at the left boundary. Now the holes 
are injected at rate 7 and removed at rate a at the left boundary and 
they are injected at rate /? and removed at rate 5 at the right bound- 
ary. They also jump with the same exclusion rules as the particles in 
the bulk. Therefore, this symmetry implies that 

^x{a, 7, 5, P, z, N) = /x(7, a, P,d,^,N) . (2.8) 

3. The Gallavotti-Cohen symmetry: the rates a, (5, 7, 5 represent the trans- 
fer of particles between the system and reservoirs at densities pa = 

and = at the two boundaries (sites 1 and N). When pa = Pb, 
the system is in equilibrium and the dynamics satisfy detailed balance 
with respect to a Bernoulli measure [20] at density p = pa = Pb- One 
can always think that the case pa / Pb represents the effect of an ex- 
ternal field which enhances the flux of particles from one reservoir into 
the system, a situation for which (as explained in the Appendix A) 
the Gallavotti-Cohen relation holds [32, 33]. This implies that 

fi{a, 7, S, 13, z, N) = p{a, 7, S,(3,^,N). (2.9) 

apz 

2.3 Main results 

When N is large, one finds, at least pertubatively in powers of a and 7 (see 
sections 5 and 6) that p, depends only on the densities pa and pb of the left 
and right reservoirs 

Pa = ^— ; Pb=lT-T 2.10 

a -I- 7 p + d 

instead of the four parameters a, (3, 7 and 5. 

The three symmetries (2.7)-(2.9) then become 

Kpa,Pb,Z,N) = p{pb,Pa,K^) (2.11) 

p{pa, Pb, z, N) = p{l - p„, 1 - Pb, Kn) (2.12) 
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KPa, Pb, Z,N) = l^i Pa, Pb, -^^^4 



n v-y (2.13) 

^Po(l - Pb) J 

It is also a fact observed in the perturbation theory to arbitrary order (see 
sections 5 and 6) that, for large N, pis proportional to 1/iV times a function 
of a single variable lo defined by 



CO 



{Z - ^){PaZ - Pb- PaPb{z - 1)) 



(2.14) 



and the result of our expansion in powers of lo of section 6 is that 

^ = lHM + o(i,; 



(2.15) 

where 

The symmetries (2.11)-(2.13) leave uj given by (2.14) unchanged so that p 
given by (2. 14), (2. 15) satisfies automatically these symmetries. From (2.16) 
one can easily obtain the large A'^ expression of the first four cumulants of 

Q{t) 



lim 

t— >oo 


m)) ^ 
t 


^[Pa 


- Pb] 






lim 

t— >oo 


t 


1 

~ N 


Pa + Pb- 


'^{pI + PaPb + Pb) 

3 




lim 

t—*oo 


{QHt))c 
t 


- ^{Pa-Pb) 


15 


lim 


{QHt))c 
t 


1 

~ N 


Pa + Pb- 


2{7pl + PaPb + 7pl) 
3 



(2.17) 
(2.18) 

(2.19) 

(2.20) 



32p^ + Spipb + 8paPt + 32pl 96pl + Gip^pb - 40plpl + ^"^Papt + 96p^ 
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One can notice that the nth cumulant (at least for n < 4) is a polynomial 
of degree n in Pa,Pb- We will comment on this at the end of section 5. 
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2.4 Two particular cases 

Let us examine these expressions in two particular cases. First when 

Pa = l and pb = (2.21) 
one sees that (2.17)-(2.20) give in the long time limit 

{Q{t)) , 1 



t N 

{Qmc , 1 

t ^ 3A^ 
{Q{t?)c , 1 
t ^ 15N 

{Qit)% , -1 



(2.22) 
(2.23) 
(2.24) 
(2.25) 



These numbers coincide with those found for quantum conductors with many 
channels in the metallic regime [5] and for quasi-classical conductors anal- 
ysed by a Boltzmann-Langevin approach [10]. 

Another particular case of interest is when the two reservoirs are at the 
same density p 

Pa = Pb = P ■ (2.26) 
All the odd cumulants vanish and (2. 18), (2. 20) become 

(Q^ ^ 1 _ ^2.27) 
Inn M^lk ^ i.2p(l - p)(l - 2pf (2.28) 
2.5 Conjecture 

We see in (2.28) that the fourth cumulant vanishes when Pa = Pb = 1/2- 
We conjecture that in this particular case, Pa = Pb = 1/2, all the higher 
cumulants vanish (i.e. the distribution is Gaussian) so that p is given in 
this case by 

l(log^ 
^ N 4 ^^\N^ 
This conjecture (see (2. 14), (2. 15)) fully determines the function R[u)) 

R{u) = [log (ym^ + V^)] ^ (2.29) 
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and therefore /x, using (2. 15), (2. 29) for arbitrary Pa,Ph and z. 

Expression (2.29) needs to be modified when uj becomes negative. We 
will also conjecture that for a; < 0, (2.29) is replaced by its analytic contin- 
uation 

R{uj) = - [sin-^ {y-^y\ ■ (2-30) 

Looking again at the first case we analyzed (2.21), we get that not only 
the first four cumulants (2.22)-(2.25) are the same as those of the distribution 
first obtained by Lee, Levitov and Yakovets [5] , but all the higher cumulants 
are the same 

{Qmc , -1 

t ~^ 105N 

{Qmc . 1 

t 231N 
{Q{ty)c 27 



t 5005iV 

{QitfU , -3 

t ^ 715N ■ 

In the equilibrium case {pa = Pb = p) too, this conjecture determines all 
the cumulants higher than (2.27), (2.28) 

^ 2p(l - p){2p - 1)2(1 - IQp + 16/) 

^'^^^^^^^ 2p{l - p){2p - 1)2(1 - SOp + 656/32 _ ii^2p^ + 575^4) _ 

Our conjecture for the distribution of Q{t) for arbitrary pa and ph co- 
incides with the distribution found in a multi-channel quantum picture [7] 
(with a small discrepancy with the distribution proposed in [6]). 

2.6 The large deviation function 

The knowledge of the large behaviour of p gives some information on 
the large deviation function FN{q). This large deviation function -F/v(q') is 
defined by 

Probability ~ ~ exp [tF^(g)] (2.31) 

or for a more mathematical definition 

lim - log [Probability {tq < Q{t) <tq + 1)] = Fjv(g) (2.32) 

t—*oo t 
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If we knew fi{a, (3,5,'y, z, N), one would determine Fj^{q) in a parametric 
form by varying z 

djj, 



oz 



FN{q) = fi- qlogz . 

As here we know only jj, only for large A'' (see (2.15)), we cannot get the 
full large deviation function -Fiv(g) for arbitrary N but we can say that in 
the large N limit, 



where the function G{q) can be constructed from R{u) in a parametric form 
by varying z 

' duj\ f dR{ijj) 



dz J \ duj J ^ ^ 

G{q) = R{u) - q log z (2.35) 




Figure 1: The rescalcd large deviation functions G{q) versus q in the cases 
Pa = Pb = 1/4 (left thick curve) and pa = 1 and pi, = (right thick curve). 
The thin lines represent for comparison the Gaussians with the same two 
moments. 
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This means that for large N we know -F/v(g) only for deviations q of order 
1/N. Figure 1 shows G{q) versus q for two choices of pa and pb (the case 
Pa = 1 and Pb = and the case Pa = Pb = -25). 

3 The average current and its variance 

In this section we show that the expected value and the variance of the 
integrated current Q{t) can be calculated easily by using the conservation 
rules. 

Let us define Yi{t) the integrated current between sites i and i + 1 during 
the time interval 0, t (so Yi {t) is the total number of particles which have 
jumped from i to i + 1 minus the number of particles which have jumped 
from i + 1 to i during time t). Similarly let us define lo(^) the integrated 
current from the left reservoir to site 1 and Yj\f{t) the integrated current 
from site N to the right reservoir. Note that Yo{t) and Q{t) have exactly 
the same definition and therefore 

Q{t) = Yo{t) . 

The conservation of the number of particles implies that 

Yi{t) = Yi_i{t) + Ti{0)-Ti{t) . (3.1) 

The difference between Yi{t) and Yj{t) remains bounded ((3.1) implies that 
\Yiit) - K,_i(t)| < 2 and \Y,{t) - Y^{t)\ < 2\j - i\ ). Therefore in the long 
time limit the cumulants of Yi{t) do not depend on i. 

Um '.^^ = lim "'^<-''""> = li,n . (3.2) 

t-^OO t t—^00 t t— >00 t 

The very definition of the dynamics in section 2 means that during each 
time interval dt <^ 1, 

Yo{t + dt) = Yo{t) with probability 1 — a{l — Ti)dt — ^T\dt 
Yoit) + 1 with probability a{l — Ti)dt 
Yfjit) — 1 with probability ^T\dt 

Prom this evolution one can deduce the following time evolution for the 
moments of YQ(t): 

— - — = a-(a + 7)(ri) (3.3) 
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^^^^^ = 2a{Yo{t)) - 2{a + i){Y^{t)T{) + a + (7 - a)(ri) (3.4) 
and more generally 

^^^^ = a([(FoW + l)'^ - Yo{tf]{l - n)) + 7([(lo(t) - 1)'= - Y^{tf]T,) 
Prom (3.3), (3.4) we obtain 

j^{yo{tf) - {y^{t)? = -2(a + 7)[(io(t)ri) - (yo(i))(ri)] + « + (7 - a)(ri) 

(3.5) 

Similarly starting from the dynamics of the integrated current Yi{t) 
through the bond i,i + 1 or of the integrated current Ypf{t) between site 
N and the right reservoir one can get 

— — = in) - {n+i) (3.6) 

^(y,(i)')-(li(t))' = 2[(y,(0(ri-T,+i))-(yi(t))(ri-ri+i)]+(ri+T,+i-2riri+i) 

(3.7) 

and 

^JM)1 = ^^ + S)ir,)-S (3.8) 

|(yAr(t)2) - (yiv(t))' = 2(/3 + 5)[(yAr(t)rAr) - (Fiv(i))(rAr)] + <5+ (/3-<5)(rAr) 

(3.9) 

3.1 The current 

If we define the parameters a,b,pa,Pb as in [19, 20] and (2.10) 

1 . 1 



a+7 p+S 

a 5 



(3.10) 



Pa = \ ; Pb 



a+7 P+S 
one obtains by combining (3. 3), (3. 6), (3. 8) that 



d{Yo) ^MYn) ^^'d(y,) 

1=1 
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We know (3.1,3.2) that in the steady state d{Yi)/dt does not depend on i. 
Therefore, one obtains that way the steady state current 

d{Q) ^ djyj) ^ Pa-Pb 

dt dt N + a + h-1 ^ ' ' 

which gives (2.17) for large TV". 

3.2 The variance 

Similarly adding (3.5), (3. 7), (3. 9) and using the fact that in the steady state 
\t^'^ ^'^^^ i^ot depend on i one gets 

{a + h + N-l) = 2Y.tl{y^r^) - {Y,){n) - {Y^.^n) + (Fi-i)(Ti) 

(3.12) 

and using (3.1) one obtains (using that (rj(O)ri(t)) (ri(0))(ri) in the long 
time limit) 

i=i i=l 

(3.13) 

All the steady state correlations can be calculated exactly [34, 20], in par- 
ticular 

, , , N-i + h . . pa{N + b-i)+ pb{i-l + a) 
^^^^ = P'^ N + a + b-& '^ - = N + a + b-l 

and for i < j 

,2 {a + i-l){b + N-j) 



inTj) - {Ti){Tj) = -{pb - paY 



{N + a + b- l)2(Ar + a + 6 - 2) 
so that (3.13) becomes 

dm - {Q}'] d[{Y^) - (y,)2] 1 

It = It = + - 

a(a - l)(2a - 1) + 6(6 - 1)(26 - 1) - Ni{Ni - l){2Ni - 1) , 
+ 37V3(iVi - 1) " 

(3.14) 

where Ni = N + a + b— 1. In the large N limit, one obtains (2.18). 
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4 A hierarchy of equations for the correlation func- 
tions 

In the long time limit, the vector Vt{C,z) in (2.4) becomes an eigenvector 
of the matrix Mz defined in (2.5). 

where i^^iC) satisfies 

/xV/.(C) = ^M,(C,C')Vm(C') (4.1) 

c 

From (4.1), one can build a hierarchy of equations which, as we shall see 
it in the next section 5, can be truncated either when one expands in powers 
of 2; — 1 to obtain the first cumulants or when the densities pa and in the 
reservoirs are small. 

Let us define the following correlation functions: 
for 1 < i < iV 

Ti = Y,i^^{C)n{C) (4.2) 
c 

ioi l<i <j <N 

Ui,j = Y.MC)n{C)Tj{C) (4.3) 
c 



foTl<i<j<k<N 



Vi,j,k = Y^MCMC)Tj{C)Tk{C) (4.4) 
c 

and so on, with the convention that 

c 

Inserting these definitions into (4.1), one obtains a hierarchy of equations 
for the one-point functions Tj, the two point functions Uij, and so on. By 
summing (4.1) over all C, one obtains 

fi = a{z — 1) + ^7 az + a — ■y^ Ti (4.6) 
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By multiplying (4.1) by Ti{C) and summing over C, one gets 

IJ,Ti = a{z - l)Ti + (^j^-az + a-j^ Ui^i + Ti_i + T^+i - 2Ti (4.7) 

At the two boundaries (4.7) is modified 

/xTi =az- {az + -f)T^ + T2-T1 (4.8) 

liTN = a{z - 1)Tn + (^7^ - + a - 7^ Ui,n + Tjv-i - (1 + /3 + 5)rAr + 6 

(4.9) 

In fact (4.8) and (4.9) (which are the boundary versions of (4.7)) reduce 
to (4.7) provided that we require that Tq, T/v+i and Ui^i (for non-physical 
values of the parameters) satisfy 

a{z - l)Ti + ^7^ -az + a-j^ J7i,i + To-Ti = az- {az + 'y)Ti (4.10) 

S-{I3 + 5)Tn = Tn+1 - Tn (4.11) 
Similarly by multiplying by Ti{C)Tj{C) one gets 

nUi^j = a{z - l)Uij + ^7^ -az + a-jj Vi^ij + Ui-ij + C/j+ij 

+ - 4f/i,,- (4.12) 

the boundary conditions and the case j = i+1 being automatically satisfied 
provided that the extensions of J7o,i) Ui^i, Ui^N+i, to non-physical values 
satisfy 

a{z - l)Ui^i + (^1^ -az + a-'jj -I- Uo,i - Ui^i = azTi - {az + 7)C/i,i 

(4.13) 

6Ti -{P + 6) Ui,N = Ui,N+i - Ui,N (4.14) 

Ui^i + Ui+u+i = 2Ui,i+i (4.15) 

(Note that definitions such as (4.3) do not tell us what Ui_i is as Uij is only 
defined for j > i. In general the values one has to choose for non-physical 
values of the parameters in order to satisfy the boundary conditions (4.13)- 
(4.15) are different from what one could obtain by simply putting j = i 
in the definition: in particular Ui^i ^ T^. In this whole paper the C/ij we 
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calculate are always polynomials in i and j and the unphysical values such 
as Ui^i are simply obtained by taking j = i in the polynomial Uij). 

All the relations for the higher correlation functions can be generated 
along the same steps. One way of writing all these relations is to introduce 
the generating function 



$(ai,...aL;z) 



z'^Wexp 



\aiTi{t) 



(where as above, is the total number of particles transferred from the 

left reservoir to site 1 during time t) . For large t one expects that 



$(ai,...aL) ~ e''* ^(ai,...az,) 



where 6 satisfies 



n4) 



L-l 



^2 



i=l 



+(e"'-"'+i - 1) 



a 92 



dai+i dttidai+i J 



1 



dai 



+S (e"^ - 1) 1 - 



d 



+ p (e-"^ - 1) 



d 
daL 



(4.16) 



Expanding (4.16) in powers of the allows one to recover all the above 
relations between the correlation functions (4. 6), (4. 7) and to generate 
the equations satisfied by the higher correlations. The first levels of the 
hierarchy are summarized in Appendix C. 



5 The low density expansion 

When the densities pa and pi, of the reservoirs are small the n-point function 
is of order n to leading order in pa and p^. To calculate /x to order n in pa and 
Ph, one can truncate the hierarchy by neglecting all the m-point correlation 
functions for m > n. 

A priori the truncated hierarchy remains a problem hard to solve. How- 
ever we noticed that the solutions Ti, Uij ... of the truncated hierarchy are 
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always polynomials in the coordinates i,j.. (see the Appendix B on the anal- 
ogy with a multi-particle ruin problem). For example if one tries to expand 
to order 3 in pa and pi, (or in a and 7) one finds that Tj is a polynomial of 
degree 5 in i, Uij of degree 4 in i,j and Vi^j^k a polynomial of degree 3 in 
i,j, k (in fact Vij^k is linear in each of the three coordinates). So to solve the 
truncated hierarchy, we introduced arbitrary parameters (the coefficients of 
all the polynomials in i, i,j, i,j,k ...) and the equations of the hierarchy 
give us a finite set of linear equations to solve for these parameters. 

We used Mathcmatica to solve these linear equations. The expressions 
become quickly complicated for general a,b,pa,Pb- The general expression 
(C.14) of /z at order 3 in and pi, is given in the Appendix C. We give here 
the result obtained that way for p to order 3 in pa, pi, when a = b = 1 



(5.1) 



(Z - l){paZ - Pb) 

z{N + l) 

{Z - I fjpl + ApaPbZ + pIz^ + 2N{pI + PaPbZ + p^Z^)) 
6z2(iV + l)2 

{z - lf{2N + l){paz - pb){3pl + 9paPbZ + Splz^ + iV(4pg + 7papbz + iplz^)) 

45^3 (iV + 1)3 



For large N the expression (C.14) of p gets much simpler: to leading 
order in N, the results do not depend anymore on the two parameters a and 
b and one gets 



N 



{pgZ - pb){z - 1) {p^Z + PaPbZ + Pf,){z - 1) 



3^2 



2(p„Z - Pb)Mz'^ + 7paPbZ + pI){z - 1)3 

45^3 



(5.2) 



So in this large N regime, p is proportional to 1/N and is a function of three 

parameters Pa, Pb and z. In fact, if one uses the parameter uj = {z — \){paZ — 
pb — papb{z — 1))/ z defined in (2.14), one can easily check that (5.2) can be 
rewritten as 

^c.2 + -|^=^ + 0(c.^)) (5.3) 



Up to the factor 1/iV, p depends on the single parameter uj, defined by (2.14), 
(at least to order 3 in u;). The expansion of p in powers of pa and pb to 
third order determines exactly the first three cumulants and more generally 
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the expansion of /j, to order n would give the exact expression of the first n 
cumulants. This can be understood by noticing the similarity between an 
expansion of n in powers of z — 1 and an expansion of /i in powers of pa and 
pi,. In both cases, the hierarchy can be truncated and one can neglect all 
the correlations higher than the ra-point function if one wishes to obtain p 
at order n. This is the reason why the exact expression of the nth cumulant 
is a polynomial of degree n in pa and pb as noticed at the end of section 2.3. 

6 Continuous limit 

The expressions of the Tj's, C/jj's , VJj^fc's we have obtained by Mathematica 
to solve the hierarchy in powers of pa and pb are rather complicated. However 
they take a somewhat simpler form in the large N limit. If one considers 
the connected correlation functions Uij,Vij^k ■■■ defined by (C.15), (C.16), 
their expressions become functions of the continuous variables: 

To leading order in and to third order in powers of pa and pb, one 
obtains that way: 

Ti ~ paZ [rxi + (1 - xi) (l + s(-l + (1 - rfxi/3 - (1 - rfxf/G) 
+s^{l + (-23 + 24r - Qr^ + 8r^)xi/45 + (29 - 42r + 27r'^ - 14r^)x?/90 

+(r - lfxl/15 - (r - 1)^x^/60)) 1(6.2) 



Uij ~ ^pIz'^ [xi{l - X2) (-(1 - rf + s(2(4 - 6r + Sr"^ - r^)/S+ 

(r - 1)^X1 - (r - 1)^x^/3 + 2(r - 1)^X2/3 - (r - lfxl/3))] (6.3) 



^i,J,k ^ -^py [-2xi(l - 2x2)(l - xs)] (6.4) 



where the parameters r and s are defined by 



r = (6.5) 

PaZ 

S = pa{z - 1) . (6.6) 
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By examining the hierarchy (see Appendix C) for the connected functions 
and by assuming that the structure obtained up to third order persists to 
higher orders one expects that to leading order in N 

Ti ~ pazfix,) (6.7) 

2 2 

Ui,j^^g{xi,X2) (6.8) 

■^M.fe - -^h{xi,X2,X3) . (6.9) 

With this scahng and as /it is of order 1/N, one can neglect the right hand 
side of (C.23)-(C.25). One can even show that g{0,X2) = h{0,X2,X3) = 
i{0,X2,xs,X4) = so that (see (C.22)) 

Ti ~ (6.10) 

az + J 

e Ui^i ~ ^1 - (1 + Pa{z - l))('Ul,i - UQ,i) (6.11) 

e zi^ij c=^(l-^^{l + pa{z - l)){zi,i,j - zo,i,j) (6.12) 

and with these simplifications the hierarchy (C.18)-(C.31) in the large N 
regime becomes: 
the equation for /x 

p = sil + s)f'{0) (6.13) 



the bulk equations (C.19-C.21) 



S(l + S) — /l(0,X2,X3) = ( ^ ^ ) 3{X2,X3) (6.15) 



d ( d'^ d'^ d"^ 

S(l + s) — i{0,X2,X3,X4) = ( ^ + ^ + ^ ) h{x2,X3,X4) (6.16) 

the left and right boundary equations 



/w = TT^ ; /w = ^ (6.17) 
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g{0,x)=g{x,l)=0 
h{^,x,y) = h{x,y,l) = 
the equations for adjacent particles 



d 
dxi 



d 



g{x,x) 



dX2, 

h{x,x,y) 



_d d_ 

dxi dx2 , 

4— - h{x,y,y) 

dx2 dxz J 



-2 



-2 



df{x) 
dx\ 

dfjx) dg{x, y) 
dxi dxi 

df{y)dg{x,y) 
dxi dx2 



d 


d 


dxi 


dX2 


d 


d 


dX2 


dX3 


d 


d 


dX3 


dX4 



i{x,x,y,z) 



^ df{x) dh{x,y,z) ^dg{x,y) dg{x,z) 



1 1— I iix,y,y,z) 



dxi dxi 
, df{y) dh{x,y,z) 
dxi dx2 
df{z) dh{x,y,z) 



- 2 



dxi dxi 
, dg{x,y) dg{y,z) 

dx2 dxi 
dg{x, z) dg{y, z) 



(6.18) 
(6.19) 

(6.20) 

(6.21) 
(6.22) 
(6.23) 
(6.24) 
(6.25) 



dx\ dxs dx2 dx2 

One can then solve this hierarchy, up to an arbitrary order in s. To find 

H at order n in s, one needs to know / at order n — 1, g at order n — 2 and 
so on. When we calculated /x at order 4 in s, we obtained 



1 



(xi - 1) 



1 



xi{xi-2) 2 (3a;? -12x1 + 28^1 - 32) 
to h UJ 



6 180 

^xi{5xl - 30xj + 138x1 - 352a;? + 600a;i - 576) 



5040 



(6.26) 



g{xi,x2) 



Xi{x2 - 1) 



l-U 



2 - Sxi + xf - 2x2 + x\ 



X^ 2,3 
— - + 

24 4 



xl{26 - 10x2 + 5x1) a;i(3- 2^2 + a;i) 56 - 80a;2 + 60a;| - 20a;| + 5a;^ 



36 



+ 



120 



(6.27) 



h{xi,X2,X3) = r 



s + 1 



a;i(a;3-l)[4a;2-2+ 



—iO 



5^2 



lbx\ + 5x2 (xf — 3x1 + a;3 — 2x3 + 4) — 3x\ 



9x1 



2x3 



4x3-6 



6.28) 
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i{xi,X2,X3,Xi) = (r- 



s + 1 



2xi{x4 - 1) (15x22:3 - 10x2 - + 3) (6.29) 



and 



1^ 



N 



a;2 8cj3 4cj4 
— + — — + (u^) 



3 45 35 



(6.30) 



where we give the expressions of f,g,h,i and fi (except for a simple factor 
r — j^) in powers of uo defined by (2.14) instead of s. 

In principle all the expressions should depend on the two parameters 
s and r, but we observe that they only depend on the single parameter 
CO. This can be understood by noticing that if f,g,h... solve the hierarchy 
(6.13)-(6.25) for a certain choice of r, s, then Af + B, A^g, A^h ... solve the 
same hierarchy for r', s' with the same value of fi if r' and s' satisfy 



1 + s' 



A- 



l + s 



+ B 



r' = Ar + B 
As'{l + s) = s{l + s) 
These three relations are compatible only when 



r's' 



rs — rs 



so that when uj = s — rs — rs^ remains unchanged, one can easily transform 
the solution of the hierarchy leaving /x unchanged. 



7 Conclusion 

In this paper we have obtained the first four cumulants (2.17)-(2.20) of the 
integrated current for the symmetric simple exclusion process with open 
boundaries. To our surprise, the generating function of the integrated cur- 
rent (2. 6), (2. 15) depends on the densities of the reservoirs pa and and on 
the fugacity z, the parameter conjugated to the integrated current, through 
a single parameter oj defined in (2.14). It would be interesting to understand 
why this is so through a simple physical argument. 

When /9a = /95 = 1/2, the fourth cumulant vanishes and we have conjec- 
tured that in this particular case^ the distribution of the integrated current 
Q{t) is Gaussian (in the range ~ j^). Based on this conjecture, we can 
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predict (2. 29), (2. 30) the large deviation function of the current for arbitrary 
choices of Pa,Pb- For Pa = 1 and pt = 0, the distribution of the integrated 
current we obtained is identical to the one known for one dimensional quan- 
tum conductors in their metallic regime [5, 4]. 

The similarity between these results is striking if we consider the drastic 
differences in the corresponding formalisms. In the quantum treatment of a 
diffusive conductor, the statistics of the time integrated current appears as 
the result of a convolution of a large number of independent binomial laws, 
one for each conduction channel [5]. In the limit of a large number of such 
channels (i.e. when the transverse dimension of the conductor is much larger 
than the Fermi wave length) the result of this convolution is governed by the 
universal distribution of eigenvalues of the transmission matrix for a single 
particle in the presence of quenched disorder. The exclusion effects induced 
by the Pauli principle only appear in the selection of the energy window 
in which single particle states contribute to the current. By contrast, the 
classical model considered here has no transverse degree of freedom, and 
the exclusion constraint plays a crucial role. To our knowledge, a complete 
understanding of the connection between the two models is still lacking. We 
simply conjecture that an intermediate description in terms of a Boltzmann 
equation with additional noise terms, as developed for instance in [8, 10] 
for the quantum diffusive case, may help to bridge the gap between the two 
classes of systems. 

The first open question left at the end of the present paper is whether one 
could prove or disprove our conjecture for p in section 2.5. It would also be 
interesting to sec the degree of universality of the results obtained here, i.e. 
how much they depend on the precise definition of the model. In particular 
it would be nice to see whether a more macroscopic approach could be used 
to calculate the fluctuations of the current [21, 22]. Another open question 
would be to know how our results would be modified by an asymmetry 
[18, 36, 37] in the bulk, in particular in the case of a weak asymmetry [38]. 
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A The Gallavotti- Cohen relation 



In this appendix we rederive, following Lebowitz and Spohn [33], the Gallavotti- 
Cohen relation for a system with stochastic dynamics in contact with several 
reservoirs of particles. 

Let us consider an irreducible Markov process for a system with a finite 
number of internal configurations C. We assume that this system is in con- 
tact with a reservoir A (or several reservoirs A,B,C) and that during each 
infinitesimal time interval dt, there is a probability Wq{C' ,C)dt of a jump 
from C to C with q particles transferred from reservoir A to the system 
during this jump. As the system is in general in contact with other reser- 
voirs, these particles might later on be transferred to other reservoirs, so 
that Wo(C',C) allows jumps where the number of particles in the system is 
not conserved. 

Imagine that the system is in equilibrium with reservoir A, that is the 
jumping rates Wq{C',C) satisfy the detailed balance condition 

WgiC',C)P,^iC) = W^_,(C,C')Peq(C') (A.l) 

where Peq(C) is the steady state probability of the Markov process. 

Clearly the detailed balance condition (A.l) implies that the average 
current of particles vanishes and that the probability of seeing any given 
jump is equal to the probability of its time reversal as it should for a system 
at equilibrium. 

Now let us modify the dynamics by introducing a field E which enhances 
the injection of particles into the system so that Wq{C',C) is replaced by 

e^'^Wg{C',C) (A.2) 

This field E produces a current which of course fluctuates due to the stochas- 
tic nature of the Markov process. 

Let us denote by Q{t) the total number of particles transferred from 
reservoir A to the system during time t and Rt{C,Q) the probability of 
Q{t), given that the system is in configuration C at time t. The evolution of 
Rt{C, Q) is clearly 

^Rt{C, Q) = E E [Wg{C,C')Rt{C', Q-q)- WgiC C)Rt{C, Q)] 
If one introduces the generating functions rj (C, A) defined by 

n(c,A) = ^e^«i?t(c,g) 

Q 
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they evolve according to 

^MC, A) = 5^5^ {e(^+^)W,(C,C')rt(C', A) - e^W,(C',C)n(C, A)} 
q C 

This impUes that for large t, 

where iJ,{X, E) is the largest eigenvalue of the matrix M\ e 

Mx,E = ^e(^+^)'' W,{C,C') - 5(C,C')EEe''' W,{e",e) (A.3) 



where = 1 if C = C and if C / C . Therefore to obtain /x(A, E), one 

has to find either the right eigenvector '^r{C) of this matrix which satisfies 

/x(A, E)i^R{C) =J2I1 e^""^^^' WgiC,C')i^RiC') - E E ^li^'^ C)V'i?(C) 

(A.4) 

or its left eigenvector ipiiC) 

q C q C 

(A.5) 

Now if we use the detailed balance condition (A.l) for the first term in 
the r.h.s. of (A.5), we get 

q C ^eql,L-j ^ ^, 

(A.6) 

This shows that V^£,(C)Peq(C) is the n^/ii eigenvector of the matrix M_x-2E,e 
defined in (A.3). 

So the matrices M\^e and M-\-2E,e have exactly the same eigenvalues. 
In particular this shows that 

lx{\,E) = n{-\-2E,E) (A.7) 

which is the Gallavotti-Cohen relation. 

In the symmetric exclusion process, as described in section 2.1, we know 
[20] that detailed balance is satisfied whenever 

oi 5 / A o\ 
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If we fix P and 5 and vary a and 7, the detailed balance condition (A. 8) 
is no longer verified. However, one can always think of the variation of a 
and 7 as the effect of an external field E trying to enhance the number of 
particles transferred from the left reservoir to site 1. If one writes 

a = a'e^ and 7 = j'e"^ 

with 
and 

7(5 ' 

one sees that the system satisfies detailed balance for a',^',fj,6. Therefore 
the Gallavotti-Cohen symmetry implies (2.9) for the SSEP. 
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B The analogy with a multi-particle ruin problem 



In this appendix, wc show the similarity between the equations one has to 
solve at each level of the hierarchy and the equations which one can write 
in a multi-particle ruin problem. 

Consider first a single particle which diffuses on a chain of N sites with 
open boundary conditions. If the particle is at site i at time t, it jumps, 
during an infinitesimal time interval dt, to site i + 1 with probability dt (for 
1 < i < A^— 1) and to site i—1 with probability dt (for 2 < i < N). Moreover, 
a particle at site 1 is absorbed at the left boundary with probability adt and 
a particle at site A'^ is absorbed at the right boundary with probability f3dt. 
In the usual ruin problem [39], one asks the following question: what is 
the probability that a particle starting at site i will escape at the left 
boundary. Clearly Tj satisfies for 2 < i < N — 1 

Ti+i + Ti_i -2Ti = 

and at the boundaries 



a + T2-(l + a)ri = 

Tn-1 - (1 + f3)TN = 

These are precisely the equations (4.7)-(4.9) we had to solve in section 4, if 
one takes 'j = 6 = and z = 1 (which implies that = see (4.6)). 
The solution of this ruin problem is of course linear in i 

Ti = . (B.l) 

Let us now generalize the ruin problem to two particles (the generaliza- 
tion to more particles is straightforward). Consider two particles initially 
at sites i < j which diffuse in the same way as in the one-particle ruin 
problem, except that the two particles are not allowed to occupy the same 
site. As time goes on, one of the two particles will escape at one of the two 
boundaries, then the other particle will diffuse until it also escapes. 

Now we want to calculate the probability Uij that both particles will es- 
cape through the left boundary. One can write down the equations satisfied 
by Uij 
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Ui_ij + + + Uij+i - 4[/,,j = (B.2) 

Ui^i + f/i+i,i+i = 2Ui^i-^-l (B.3) 

(3 + a)Uu = aTi + U2,^ + + (B.4) 

(3 + (3)Ui,N = Ui^N-l + f/i+l,7V + f/i-l,7V (B.5) 

and they are identical to (4.12)-(4.15) when 7 = 5 = and z = 1 (implying 
that n = 0). It is not obvious a priori that the solution of these equations 
is simple. However the solution turns out to be linear in i and j 

(iV + l + i-l)(iV + i + i-2) • ^^-'^ 

We also see in (B.1),(B.6) that the correlation between the two particles is 
weak ^ ^ 

= Ui^j - TjTj = i — 



N 



when i and j are of order N. This weak correlation [34], which are similar 
to those seen in (6.8), (6.9), is however responsible for the non-Poissonian 
character of the fluctuations of the integrated current. 

Another quantity which has a simple expression (i.e. for which the so- 
lution is linear in i and j) is the probability Ui^j that one particle escapes 
at the right and the other particle escapes at the left, without specifying on 
which side the first particle to escape leaves (starting with two particles at 
positions i and j). In this case the equations to solve are again (B.2), (B.3) 
with boundary conditions (B.4), (B.5) replaced by 

(3 + a)Uu = a(l - Ti) + U2,^ + C/i,»+i + (B.7) 

(3 + p)Ui,N = m + Ui,N-i + Ui+i,N + Ui-i,N (B.8) 
and the solution is 

_ A^(/ + j) 2ij 2N + 2/ + 1(2A- - j - f) + l^jj + j - 2) + ^ 
(^-l + ^ + ^)(^-2 + ^ + ^) 

If one asks however a slightly more precise question, namely what is 
the probability Uij that (starting with a particle at i and a particle at 
j), the first particle to escape leaves at the right boundary, and then the 
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remaining particle escapes at the left boundary, the equations to solve are 
still (B.2),(B.3) but with the boundary conditions (B.4),(B.5) replaced by 



(3 + a)Ui,i = U2,i + Ui,i+i + ?7i,i_i 



(B.9) 



(B.IO) 



These new boundary conditions make the problem much harder and one can 

check that the solution is no longer linear (or even polynomial) in i and j. 

So the same problem (B.2),(B.3) with the boundary conditions (B.4),(B.5) 
or (B.7),(B.8) is easy (the solution is linear in i and j) whereas it is hard 
with boundary conditions (B.9), (B.IO). The main reason which made pos- 
sible the calculation of the cumulants in the present paper is that each time 
we had to solve equations of the type (B.2),(B.3), the boundary conditions 
were such that the solution was polynomial in the coordinates i and j. 
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C The hierarchy 

C.l The hierarchy for the correlation functions 

The hierarchy of section 4 can be summarized as follows: 



The equation for fx 



^ = «(,_i)_(^i±2)(i^r, (C.l) 



The bulk equations 

//T, = aiz - l)Ti - («^ + 7K^-l) ^^ . ^ ^ r^,^^ _ (^.2) 

nUij = a{z-l)Uij-^ '-^ ^-Vi,i,j+Ui-ij+Ui+ij+Uij-i+Uij+i-4Uij 

(C.3) 

-|- 'yjiZ — 1) 

fJ'Vij,k = Oi{z - l)Vi,j,k Wi,i,j,k + Vi-ij,k + Vi+ij,k 

+Vi,j-i,k + Vij+i^k + Vi,j,k-i + Vij^k+i - 61^i (C.4) 

The left boundary equations 

a{z - l)T^ - («^ + 7K^-l) ^^^^ +To-Ti = az- {az + ^)Ti (C.5) 

a{z-l)Ui,i- + 7)(^ - ^^ Vi,i,i + Uo,i-Ui,i = azTi-{az + j)Ui,i (C.6) 

icxz ~I~t)(-s — 1) 

- l)Vi,i,j - -W^,i,i,j + VQ,i,3 - Vi,i,j = azUij - iaz+j)Vi,ij 

(C.7) 

The right boundary equations 

d-{P + d)TN = Tn+1 - Tn (C.8) 

STi -{(3 + 5)Ui,N = Ui,N+i - Ui,N (C.9) 
6Uij - (/? + 6)Vij,N = Vij,N+i - Vi,j,N (C.IO) 
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The equations for adjacent particles 

Vi,i,j + Vi+i^i+ij = 2Vi^i+ij (C.12) 

+ ^jj+ij+i = 2V^ij,j+i (C.13) 
When one solves this hierarchy up to order 3 in pa and pb, one obtains 



{PaZ - Pb)iz - 1) 2 
M- +1^ ij 



(a - + + 6 - 36^ + 2b^){paZ - pbf 

6z27Vf(iVi - 1) 

[PaZ-pbf pIz^+pI pIz^ + PaPbZ + pI 



+{z~l) 

(-7a + SOa^ - SOa^ + 4 

+ 



6z2iV2(7Vi - 1) 2z27Vi(7Vi - 1) 3z2(iVi - 1) 

(a - 3a2 + 2a^ + b- 3b'^ + 2h^f{paZ - pbf 



9z^N^{Ni - l){Ni - 2) 
(-7a + 30a2 - 50a^ + 45a* - 18a^ -7b+ 30^2 - 5053 ^ 45^4 _ i865)(p„^ - pf,f 



+ 



'i5z^Nf{Ni-l){Ni-2) 
[(2 + 15a - 45a2 + 300^ + 156 - 4562 + 3063)(p2^2 ^ ^2) _ 4p^p^z]{paZ - pb) 

Af,z^Nl{Ni - l)(JVi - 2) 
(a - 3a^ + 2a'^ + 6 - 3^2 + 2b^){plz^ - pi) + 3(p2^^ + pDjp^z - pb) 

9z3iV2(iVi - l)(iVi - 2) 

^{Plz^ - pI) jPaZ - Pb){2plz^ + ^PaPbZ + 2p2) 



+ 



Qz^Ni{Ni - l){Ni - 2) 3^3(7Vi - l){Ni - 2) 

2{paZ - Pb){^plz^ + 7papbZ + ^Pb)Nl 



+- 



A5z^{Ni - l){Ni - 2) 



(C.14) 



where Ni = N + a + b - 1, pa = a/{a + 7), pb = S/{P + 5),a = l/(a + 7) 
and 6 = + 

C.2 The hierarchy for connected correlation functions 

If one introduces the connected functions Uij, Vij^k--- defined by 

U,.j =T,Ti + n,j (C.15) 
Vi,j,k = TiTjTk + UijTk + Ui^kTj + Uj^kTi + Vij^k (C.16) 

Wij^k,l = TiTjTkTi + UijTkTi + Ui^kTjTi + uj^kTiTi + Ui^iTjTk 
+Uj,lTiTk + Uk,iTiTj + UijUk,i + Ui^kUj,i + Ui^iUj^k + ■"i,j,fc^i 

+Vi,j,iTk + Vi^k,lTj + Vj^k,lTi + w^ij.fe,/ (C.17) 
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the hierarchy becomes 



The equation for /j, (obtained by combining (C.l), (C.2) and (C.5)) 

z - 1 



■{ri - T2) 



The bulk equations 



e = 



{az + ^){z- 1) 



The left boundary equations 

{az + 'y)Ti — az + Tq — Ti = e ui^i — fi Ti 
(a + 7)ni,j + tio,j - = e vi^i^i - 2fi m^i 
(a + 7)-ui,ij + f o,i,j - = e j + iti.iMi j] - 2/1 vi^ij 
The right boundary equations 

S-{^ + 5)Tn = Tn+1 - Tm 

-{13 + 5)Ui^M = Ui^N+l - Ui^N 
-{P + S)vij^N = Vi,j,N+l - Vij^N 

The equations for adjacent particles 



(C.18) 



e ui,i = Ti_i + Ti+i - 2Ti (C.19) 
e vi^ij = Ui-ij + Ui+i J + Uij-i + Uij+i - 4uij (C.20) 

(C.21) 

where e is defined by 



(C.22) 

(C.23) 
(C.24) 
(C.25) 

(C.26) 
(C.27) 
(C.28) 



(C.29) 
(C.30) 
(C.31) 



Vi,i,j + Vi+i^i+ij - 2^1,1+1 J = -2{Ti - Ti+i){uij - Ui+ij) 
""1,3,3 + - 2vijj+i = -2{Tj - Tj+i){uij - Uij+i) 

Wi,i,j,k + '"^1+1,1+1 j,fe — '^Wi^i+i,j,k = —'^{Ti — Ti+i){vij^k - Vi+ij,k) 

-2{uij - Ui+ij){ui^k - Ui+i,k) (C.32) 



etc... 



(As already discussed right after (4.15), the (unphysical) values To,Mi.i... 
are obtained by using the explicit (polynomial) expressions of Tj, Ujj for 
i = 0,j = i). 
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